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SECTION  I 
INTRODUCTION 

Heating  a preformed  plasma  with  a laser  has  been  of  interest  because  of  basic 
physics  and  a wide  variety  of  applications.  Some  previous  calculations  have  been 
directed  toward  the  interaction  of  CO2  lasers  with  the  plasma  created  in  the  dense 
plasma  focus  (DPF)  (ref.  1).  These  calculations  showed  that  the  presence  of  small 
amounts  of  high-Z  impurities  greatly  reduce  the  electmn  temperature  achievable  in 
heating  the  DPF,  as  well  as  affecting  the  intensity  of  the  line,  recombination, 
and  bremsstrahlung  radiation  produced  by  the  plasma.  Thus,  to  treat  situations 
in  which  high-Z  materials  or  impurities  are  present,  the  treatment  of  laser  heat- 
ing and  radiation  are  highly  dependent  upon  the  level  of  ionization  in  the  plasma 
and  energy  which  goes  into  increased  ionization. 

The  laser-plasma  heating  calculations  of  McCann  (ref.  1)  were  made  with  the 
assumption  that  the  plasma  was  fully-ionized,  but  in  a separate  calculation  he 
showed  that  the  effects  due  to  ionization  are,  indeed,  quite  significant.  This 
strong  dependence  on  ionization  can  be  attributed  to  (1)  the  actual  energy  used 
to  increase  the  level  of  ionization  is  a significant  fraction  of  the  total  energy; 
(2)  most  of  the  dominant  physical  effects  in  the  plasma  are  dependent  upon  the 
level  of  ionization  in  the  plasma  or  the  average  level  of  ionization,  Z.  These 
effects  include  the  radiations  already  mentioned,  electron  thermal  conductivity, 
and  the  ion-electron  collisional  energy  exchange;  (3)  a temperature  dependent 
ionization  model  would  be  more  accurate  than  assuming  full  ionization.  One  must 
be  careful,  however,  to  use  an  appropriate  model. 

For  these  reasons,  the  corona  model  (ref.  2)  as  used  by  McCann  was  incor- 
porated into  the  calculations  of  the  laser  heating  as  performed  by  the  one- 
dimensional , Lagrangian  coordinate,  hydrodynamic  code,  DULAH.  The  incorporation 
of  both  a temperature-dependent  ionization  model  and  the  capability  to  treat  a 
variety  of  impurity  elements  has  resulted  in  a new  code,  DULAH2.  Section  II 
exhibits  the  modifications  to  DULAH  which  were  necessary  to  incorporate  these 
capabilities. 

Since  a knowledge  of  the  optical  interactions  between  the  incident  laser  light 
and  these  preformed  plasmas  would  aid  in  optimizing  the  laser  heating,  a two- 
dimensional  , optical,  ray-tracing  code,  RATRACE,  was  developed  to  examine  several 


3 


AFWL-TR-74-231 


effects  of  interest  in  these  plasmas,  namely  refraction,  reflection,  and  inte- 
grated absorption.  Some  preliminary  results  of  this  type  of  study  are  presented 
in  section  III,  and  some  one-dimensional  heating  calculations  of  exploding  wires 
are  discussed  in  section  IV. 
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SECTION  II 

CORONA  MODEL  CALCULATIONS  OF  AVERAGE  IONIZATION  LEVEL  AND  ENERGY 

The  steady-state  corona  model  is  an  atomic  model  which  assumes  that  there  is 
a balance  between  collisional  ionization  (and  excitation)  and  radiative  combina- 
tion (and  spontaneous  decay)  in  the  optically  thin  approximation,  i.e.,  all  rad- 
iative escapes  without  further  interaction.  This  model  assumes  that  changes  in 
plasma  parameters  take  place  slowly  enough  for  the  ion  populations  to  assume 
their  new  steady-state  values  at  any  instant.  Note  that  in  this  model  the  ion 
population  distribution  is  independent  of  electron  density.  The  equation  govern- 
ing these  ion  populations  can  now  be  stated  as  the  following: 


7.87x10-5  |x(Z)pji^ 


Sldil 


Where  x(Z)  is  the  total  ionization  potential  to  get  from  a neutral  atom  to  charge 
state  Z, 


L,  < I 

((Z)  -X) 


In  the  work  by  McWhirter  (ref.  2)  the  criterion  for  applicability  of  this  model 
is  given  by 

n^  < 5.6  X 103  (z+i)6  t^i/2  exp  1 1-162  x 1Q3  (Ztl.).^  (3) 

with  Tg  in  degrees  Kelvin.  Some  limiting  cases  for  this  model  for  various  values 
of  Z are  plotted  in  figure  1 derived  from  equation  (3). 

From  equation  (1)  and  the  laws  of  charge  neutrality  and  charge  conservation 
given  by 

"e  - E 'i  "i 


Figure  1.  Corona  Model  Limits  for  Various  Values  of  Z 


the  following  relationship  can  be  derived  for  the  number  density  of  species  i 

i 

He  n fi(Tj 

j=2  J ^ 


"i  = 


i=2 


(5) 


where  f.(T  ) is  eouivalent  to  n(Z+l)/n(Z)  and  is  given  by  the  inverse  of  equation 
J e 

(1). 

The  CORONA  program  (ref.  1)  does  not  properly  calculate  the  species  densities 
and  the  average  ionization  level  for  materials  with  Z <_  2.  The  explicit  forms  for 
these  cases  are  the  following  (with  the  assumption  that  a normalized  number  of 
electrons  are  present). 


no  = 1 
ni  = l/fi 

Z = fi/(fi+l) 
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(6) 
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Z = 2 


^(l+2f2] 


1 + 2f, 


_/l  + 2f2\ 
\1  + 2fJ 


In  both  cases,  the  average  ionization  energy  is  given  by 

Z 

i-1 

where  x(Z)  is  the  total  ionization  energy  required  to  reach  charge  state  Z and 
is  given  by  equation  (2). 

1.  CORRECTIONS  TO  DIFFERENCE  SCHEME 

The  equations  governing  the  ion  and  electron  internal  energies,  and  e^, 
respectively,  are  the  following: 


IT  ' ■’’i  ^ '^e1  * 32 


Tl'-PeHT-Oel 


St  - Ej 


where 


P = energy  transferred  between  internal  kinetic  energy 
dt 

Q . = energy  transferred  in  ion-electron  collisions 
e'l 

$1<  = energy  absorbed  from  laser  by  inverse  bremsstrahlung 
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E^  = thermal  heat  condition 


Q = bremsstrahl ung  radiation 


32'  + j S,  = thermonuclear  energy  deposited  in  ions  by  a's 
e 

Ej  = average  ionization  energy 


In  modifying  the  above  equations  to  include  time-dependent  ionization  effects, 
the  ionization  energy,  Ep  must  be  accounted  for  as  well  as  using  a time- 
dependent  Z in  the  expansion  of  terms.  In  some  of  the  above  terms  it  is  neces- 
sary to  derive  the  correction  terms  because  of  a nonzero  ^ Z;  these  include  the 
bremsstrahlung  radiation  terms,  the  ion-electron  collisional  exchange  term,  and 
the  heat  conduction  term.  All  of  the  other  terms  may  be  utilized  as  expounded 
by  Brueckner,  et  al . (ref.  3)  by  merely  replacing  Z with  a time-varying  T instead 
of  assuming  full  ionization. 


Using  an  implicit  difference  scheme,  the  following  are  used 


lH+i/z 


£ f-  n ^ , 

e ^e 


2aFW2  ^pi+1 


H J 


t ! 


e^n+l  _ £^n  ijn+i/2 

2Atn+i/2  Atn+1 


'tn+1/2 


Using  the  quantities  from  the  n^^  time  step  (denoted  by  n superscript)  for  the 
appropriate  zone  or  boundary  (denoted  by  half-integral  and  integral  subscripts, 
respectively),  the  appropriate  quantities  can  be  calculated  at  time  = ti^  + 
At^^l . 


The  corrections  which  must  be  made  and  the  resulting  difference  equations  are 
as  follows: 

a.  Bremsstrahlung  Radiation  Term 

The  bremsstrahlung  radiation  is  given  by 

= 03) 
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= C^'  N.  15/2  Ae^i/2 


Noting  that  N.  is  not  constant  in  time  because  the  volume  is  allowed  to  change 
although  the  total  number  of  ions  considered  is  held  constant,  and  differentiating 
equation  (14)  yields 


0 = 6 1 e . 5 Z 

1 e ^ 


Therefore,  applying  the  same  differencing  as  in  equations  (11)  and  (12), 


Q ^ • n 


1 5 . 1 

N,"  'I  Z"  'I  . " 

1 e 


The  middle  two  terms  are  terms  not  considered  in  the  previous  work  (ref.  3), 

b.  Thermal  Conduction  Term 

The  thermal  conduction  term  is  of  the  following  form: 


^ ^ I ^ 

c ~ 7 ^0  9m 


i J— --i 

IZ  in  A 9R 


8 AG  7/2] 


The  spatial  centering  of  this  equation  as  done  in  the  original  work  (ref.  3)  is 
correct.  However,  the  implicit  approximation  that  was  used  must  be  changed. 
Using  the  type  relationship  in  equation  (11)  results  in: 


This  replaces  the  approximation 

/ n+l\7/2_ 


V.  . _ 
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P T 

' I 


The  differencing  of  equation  (17)  yields. 


(■=)::: 


n+h  _ 2 ^0 


. 7 m.  1 

'S  J~'S 


a-1 


1 

n+h 

j-h 

'N'^f 

n+% 

3-h 

[(I  A„  A)j.^(Rf  - Rj:>f)4  {I  «„ 

(r;:^  - Rr“^)] 

„n4i. 

■"m 


_\7/2h+*S 

•3/2 


) 


>.[(Z  4 (Z  An  A)j.jRf‘=  . r;:^)] 

n+l\^2 


(20) 


Upon  proper  substitution  for  as  given  in  equation  (19),  this  form 

is  correct. 


2.  If’N-ELECTRON  COLLISIONAL  EXCHANGE 

The  ion-electron  energy  exchange  term  which  appears  in  equations  (9)  and  (lO) 
is  given  by  (ref.  3)  as 


N 


Qo-i  “ I - 3/t  i — " ^-i 


^0- 


fi  - ■.) 


(21) 


where 


Tg'  = I .83  X lO-io  A 


AV/2  1 


(I) 


Z2  Jin  A 


(22) 


Substituting 


"'’o'  “ ■'’o 


J7/2 


(23) 


into  equation  (21)  makes  tq"  independent  of  time  and 


Z 7/2  I 


0 


ft  ■ •-)! 


(24) 


^ t li,  ~ ^ 


ii 


i’ 

I 


* 1 


^1'^*  I~» 


Taking  the  time  derivative  of  equation  (24)  yields 


Evaluating  these  partial  derivatives  yields 


Where  the  last  term  is  a term  which  does  not  appear  in  the  original  calculations 
and  must  be  integrated  and  included  in  the  calculation  of  . 
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SECTION  III 

RAY  TRACING  STUDIES  ON  EXPLOSING  WIRES 


1.  OPTICS 

The  ray  tracing  of  a laser  beam  through  a plasma  is  accomplished  by  first 
establishing  a two-dimensional,  rectangular  coordinate  grid  which  sets  up  a 
plane  of  cells  (figure  2).  RATRACE  (ref.  4)  accomplishes  the  propagation  of  a 
single  plane-wave  ray  through  this  grid  by  the  application  of  plane-wave  optics 
at  the  cell  boundaries.  This  includes  the  law  of  refraction  at  a plane  surface 


n^.  sin  = r\2  sin  \l>2 


and  the  law  of  reflection 


(27) 


(28) 


where  n is  the  index  of  refraction  and  is  the  angle  of  incidence  with  respect 
to  the  normal;  the  subscripts  1 and  2 refer  to  the  values  of  these  quantities  on 
the  respective  sides  of  the  interface.  The  ray  is  propagated  across  the  cell  in 
a straight  line  until  the  cell  boundary  is  encountered.  Optical  effects  are  cal- 
culated at  the  cell  interface,  and  the  ray  is  propagated  across  the  appropriate 
cell.  This  process  continues  until  the  system  boundary  is  encountered.  At  this 
point  the  trace  of  the  ray  in  the  system  of  interest  is  plotted  on  microfilm. 

For  this  calculation,  the  index  of  refraction  is  defined  from  the  electron 
number  density  by 


I 


where  (cm“3)  is  the  electron  density  in  the  cell  and  (cm"^)  is  the  critical 
electron  density  which  is  defined  as  the  density  at  which  the  plasma  frequency, 

Vp,  equals  the  laser  frequency,  v^. 
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in  jji  yii]^ 


‘1'.' 
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c '■l) 


2 e 


0 nie 


(30) 


(ref.  5),  where  the  plasma  frequency  is  given  by 


(31) 


(ref.  1 ) or  by 


V = 8.98  X 10^  N„i/2(sec"i) 
P e 


(31a) 


2.  ABSORPTION 

In  the  plasmas  under  investigation,  the  path  followed  by  the  laser  light  does 
not  fully  or  adequately  describe  the  interaction.  Another  important  indicator  of 


the  interaction  is  the  amount  of  laser  light  which  is  absorbed  in  traversing  the 
plasma.  The  fraction  of  absorbed  light  is  given  by 


I 


1 - = 1 . exp 

0 


(- 1 


(32) 


where  I/I^  is  the  fraction  transmitted,  1 - I/I^  is  the  fraction  absorbed,  and 
K is  the  inverse  bremsstrahlung  absorption  coefficient  given  by 


16  I Ng2  q6  £n  A 


1 


K.  = 


a 3c  V|_2  (2m^  kT^|572  ^1  - Vp^/v^yT^ 


(33) 


(ref.  6),  or  upon  substitution  by 


1 .56  X 10-®  Z Np2  »n  A 


(33a) 
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For  use  in  RATRACE  the  absorption  irtegral  is  approximated  by  the  finite  summa- 
tion 


i for 
all 
cells 


where  K . is  the  appropriate  absorption  coefficient  for  cell  i,  and  Ar.  is  the 

a 1 1 

distance  travelled  across  the  cell. 

The  combination  of  the  optical  path  traversed  by  the  ray  and  the  amount  of 
energy  absorbed  during  this  traversal  is  used  to  characterize  the  interaction  of 
the  laser  and  the  plasma. 

3.  APPLICATION  TO  EXPLODING  WIRE  STUDIES 

Exploding  wires  produce  plasmas  which  can  be  laser  heated.  In  treating  the 
heating  of  these  plasmas,  the  inverse  bremsstrahl ung  mechanism  of  heating  is 
assumed.  The  absorption  coefficient  for  this  mechanism  is  functionally  dependent 
upon  temperature,  electron  density,  ionization  level,  and  laser  frequency.  The 
optical  effects  as  well  as  the  effects  of  the  above  parameters  can  be  simulated 
using  the  RATRACE  code. 

The  density  profiles  used  in  this  work  are  analytical  fits  to  the  data 
obtained  on  aluminum  exploding  wires  by  holographic  interferometry  (ref.  7),  where 
the  electron  density  distributions  were  found  to  be  almost  exponential ly  decreas- 
ing with  distance  from  the  center  of  the  wire. 

Experimental  studies  of  the  heating  of  exploding  wire  plasmas  are  to  be  con- 
ducted using  one  wire  and  two  parallel  wires.  Theoretical  studies  of  these  con- 
figurations have  been  made  using  the  RATRACF  code  and  have  yielded  some  interest- 
ing results  which  are  portrayed  in  figures  3 through  8.  In  these  figures  the 
contours  of  electron  density  are  all  the  same:  the  minimum  is  1.0  x 10'^  cm"^, 

the  maximum  7.8  x 10^®  cm'^,  and  they  are  spaced  1.1  x 10^®  cm"®  apart. 

Simulation  of  the  laser  beam  incident  upon  one  wire  showed  a large  amount  of 
scattered  light  and  high  absorption  for  the  CO2  laser,  X = 10.6  ym  (figure  3), 
and  essentially  no  scattered  light  and  much  less  absorption  for  the  ruby  laser. 
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X = 0.6943  M (fiaure  4).  This  is  to  be  tyoected  because  of  the  proportionality 
of  the  absorption  coefficient  to  the  square  of  the  wavelength.  The  larger 
refractive  effects  occur  for  the  CO2  case  since  the  electron  density  is  near  the 
critical  density  for  CO2  (1  x 10^^  cirr^),  whereas  in  the  ruby  case  the  electron 
density  is  far  below  the  critical  density  for  ruby  (2.33  x lO^i  cm"^). 

Simulation  of  a laser  beam  incident  upon  two  parallel  exploding  wires  used  an 
electron  density  profile  consisting  of  a superposition  of  the  two  single  wire 
profiles  separated  by  a distance  of  (2  mm).  The  same  collimated  laser  beams, 

CO2  and  ruby,  used  in  the  one-wire  case  were  applied  to  the  two-wire  simulation. 
The  results  are  shown  in  figures  5 and  6,  respectively. 

In  the  CO2  laser  case,  the  focusing  effects  of  the  two-wire  profile  are 
quite  predominant.  The  wires  function  as  a focusing  lens  with  a focal  length 
on  the  order  of  the  separation  of  the  wires  (2  mm).  With  the  same  electron 
density,  the  ruby  laser  is  transmitted  unaffected  except  for  the  4.4  percent 
absorption  of  the  beam  intensity. 

In  reality,  the  laser  beams  which  will  be  used  in  the  exploding  wire  studies 
will  not  be  collimated  but  will  be  focused  from  a 2-cm  diameter  by  a 7.5-cm 
focal  length  lens.  The  propagation  characteristics  of  this  beam  are  employed  in 
the  simulations  represented  in  figures  7 and  8 for  ruby  and  C02»  respectively. 

The  beam  is  initially  focused  on  the  point  which  bisects  the  line  joining  the  two 
wire  centers.  In  the  ruby  case  the  position  of  this  focal  point  is  unaffected. 
Note,  too,  that  this  focused  case  yields  a somewhat  reduced  absorption  as  com- 
pared with  the  corresponding  collimated  case.  This  is  because  the  beam's 
optical  path  is  through  a less  dense  region  when  focused  between  the  wires.  For 
the  CO2  laser  the  focal  point  is  moved  somewhat  nearer  to  the  lens.  This  effect 
is  expected  based  on  the  results  shown  in  figure  5 for  the  collimated  CO2  beam. 

4.  CONCLUSIONS 

Using  the  RATRACE  code  to  simulate  two-dimensional  effects  in  laser-plasma 
interactions,  some  understanding  has  been  achieved  of  the  most  effective  manner 
in  which  to  heat  a preformed  plasma  given  its  configuration  and  some  basic  assump- 
tions about  the  absorption  mechanism.  From  the  tests  which  have  been  run  to  date, 
the  CO2  laser  incident  upon  two  wires  looks  the  most  promising  in  terms  of  the 
amount  of  energy  which  can  be  deposited  in  the  plasma  by  the  beam.  From  the  pre- 
dicted absorption  and  light  paths,  some  important  points  of  comparison  with 
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diagnostic  techniques  have  been  established.  This  should  aid  In  an  understanding 
of  these  Interactions  when  these  diagnostics  can  be  performed. 
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SECTION  IV 

INITIATION  AND  LASER  HEATING  OF  EXPLODING  WIRE  PLASMAS 

In  the  consideration  of  laser  heating  of  exploding  wire  plasmas,  a number  of 

important  questions  arise  which  must  be  considered: 

(1)  In  the  initiation  phase  of  exploding  a wire  with  stored  energy  in  a 
capacitor  bank,  what  are  the  physical  properties  of  the  wire  as  a function  of 
time--phase,  temperature,  diameter,  and  density? 

(2)  What  are  the  spatial  variations  of  these  parameters? 

(3)  Given  the  shape,  power,  and  duration  of  the  laser  pulse,  what  is  the 

best  time  to  begin  heating  the  plasma?  Some  parametric  questions  which  may  also 

be  of  interest  are  the  following:  (1)  what  density  gradient  enhances  the  heat- 

ing of  a preformed  plasma  most  effectively,  and  (2)  how  does  the  laser  heating 
scale  with  pulse  length  and  energy? 

The  present  work  is  concerned  with  some  parametric,  one-dimensional  studies 
of  the  laser  heating  of  aluminum  exploding-wire  plasmas  with  a CO2  laser.  This 
study  was  conducted  in  two  phases: 

(1)  Using  the  one-dimensional  magnetohydrodynamic  code,  MAGPIE,  (ref.  8), 
consideration  was  given  to  the  initiation  of  the  wire  plasma  by  the  capacitor 
bank.  The  bank  used  in  these  studies  has  the  properties  which  are  listed  in 
table  1,  Given  these  parameters  for  the  capacitor  bank,  the  initiation  of  wires 
of  three  different  diameters  was  studied  to  determine  what  conditions  could  be 
generated  for  laser  heating, 

(2)  Given  these  initial  conditions  for  heating,  a one-dimensional  laser- 
hydrodynamics  code,  DULAH2,  was  applied  to  determine  what  heating  could  be 
expected.  This  code  was  also  used  to  look  at  the  dependence  of  heaiing,  upon 
laser  pulse  shape  and  power. 

1.  EXPLODING  WIRE  INITIATION 

Using  the  one-dimensional  magnetohydrodynamics  code,  MAGPIE,  the  initiation 
of  exploding  wires  of  various  diameters  was  studied.  These  studies  assumed  the 
parameters  for  the  stored  energy  capacitor  bank  that  are  listed  in  table  1.  The 
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Table  1 

CAPACITOR  BANK  PARAMETERS  USED 
IN  INITIATION  STUDIES 


Capacitance 

0.75  microfarads 

Resistance 

0.29  ohms 

Inductance 

1 .29  microhenries 

Voltage 

15  kilovolts 

Stored  Energy 

85  joules 

three  exploding  wires  considered  are  aluminum  wires  of  2.0-cm  length  with  diam- 
eter of  1 , 2,  and  5 mils. 

Of  obvious  importance  and  interest  in  this  type  of  calculation  is  the  equation 
of  state  used  for  the  wire.  The  equation  of  state  used  in  these  calculations  for 
aluminum  is  a three-phase  equation  of  state  (liquid,  gas,  plasma)  in  empirical 
form.  Thus,  the  sol id-to-1 iquid  phase  transition  is  ignored  for  the  purposes  of 
this  calculation.  The  solid-liquid  transition  should  occur  at  0.08  electron  volts 
with  a specific  heat  of  0.215  calorie/gram.  For  a 2.0-cm,  5-mil  diameter  wire, 
the  phase  change  requires  approximately  3.5  millijoules  of  energy.  In  light  of 
the  parameters  of  the  bank  used  in  these  studies,  this  phase  transition  is  of 
little  importance;  ignoring  it  at  this  point  is  not  significant. 

The  intent  of  this  phase  of  the  study  is  to  determine  what  are  the  best  con- 
ditions under  which  to  produce  an  exploding  wire  plasma  if  the  constraints  of  the 
capacitor  bank  are  as  already  described.  The  first  major  consideration  in  forming 
a target  for  laser  interaction  is  that  the  target,  if  possible,  be  in  the  plasma 
state.  By  examining  the  current  through  the  wire,  as  a function  of  time,  the 
state  of  the  plasma  should  be  evident  to  some  degree.  The  current  traces  shown 
in  figure  9 are  those  for  the  three  cases  which  will  be  considered.  The  dip 
which  occurs  in  the  current  is  attributed  to  the  phase  transition  from  the  liquid 
to  the  gas.  Several  important  comparison  can  be  made:  (1)  the  larger  diameter 
wires  require  more  current  to  drive  the  phase  transition;  (2)  the  larger  the 
wire  the  later  in  time  that  the  phase  transition  occurs;  (3)  the  5-mil  diameter 
wire  appears  to  be  the  largest  wire  that  this  bank  will  drive  across  the  liquid- 
gas  phase  transition;  (4)  the  current  rise,  point  of  maximum  current,  and  fall 
off  of  the  current  all  occur  on  approximately  the  same  time  scales  regardless  of 
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Figure  9 


CO  1-MIL  ALUMINUM  EXPLODING  WIRE 

A 3-MIL  ALUMINUM  EXPLODING  WIRE 

, 5-MIL  ALUMINUM  EXPLODING  WIRE 

. Exploding  Wire  Current  Traces  for  1,  3,  and  5-niil  Aluminum  Wires 


AFWL-TR-74-231 


the  wire  size  (this  is  with  the  exception  of  the  current  change  occurring  at  the 
phase  change).  Thus  it  appears  that  the  plasma  parameters  in  this  case  do  not 
significantly  modify  the  parameters  of  the  entire  circuit. 

The  purpose  of  making  these  initiation  studies  is  to  determine  the  optimum 
conditions  for  laser  heating  of  an  exploding  wire  plasma.  Some  other  indications 
of  the  state  of  the  wire  which  will  provide  useful  information  are  the  wire 
radius  as  a function  of  time  and  the  temperature  as  a function  of  time.  These 
are  shown  in  figures  10  and  11,  respectively.  Both  of  these  figures  show  some 
relationships  which  at  first  seem  strange.  There  is  obviously  no  simple  scaling 
of  these  parameters  from  one  wire  size  to  another.  Examining  figure  10  in  detail 
shows  that  the  wires  start  expanding  at  increasingly  later  times  as  the  wire 
size  increases.  This  corresponds  to  the  fact  that  the  liquid-gas  phase  change 
occurs  later  in  time  as  seen  in  the  current  traces.  Immediately  after  the  phase 
transition,  all  wires  expand  outward  with  approximately  the  same  initial  velocity 
(indicated  by  the  constant  and  equal  sloping  lines  of  the  spatial  position  versus 
time  plot).  In  the  case  of  the  3-mil  and  5-mil  wires  this  expansion  expansion 
continues  at  constant  velocity.  However,  in  the  1-mil  case,  the  expansion 
velocity.  However,  in  the  1-mil  case,  the  expansion  velocity  is  reduced  at  about 
0.7  microsecond.  At  this  point  the  magnetic  pressure  has  increased  and  is  con- 
fining the  expanding  plasma. 

In  examining  the  temperature  versus  time  plots  in  figure  11,  several  dif- 
ferent situations  exist  as  in  the  case  of  the  plasma  outer  radius.  Here  there 
is  no  simple  scaling  relationship  between  the  temperature  reached  and  the  wire 
diameter.  By  comparison  with  figure  10  it  appears  that  the  maximum  temperature 
for  each  wire  is  reached  shortly  after  the  initial  expansion  of  the  wire  begins. 
This  corresponds,  too,  with  the  transition  point  of  gas  to  plasma.  A strange 
point  is  that  the  maximum  temperature  reached  by  the  3-mil  wire  is  higher  than 
the  maximum  temperature  reached  by  either  the  1-  or  5-mil  wires.  It  appears 
that  the  1-mil  wire  has  so  little  mass  that  it  is  blown  apart  by  the  particle 
pressure  before  much  ohmic  heating  can  occur  and  raise  the  temperature.  The  3- 
mil  wire  is  held  together  longer  and  heated  more.  The  5-mil  wire  then  contains 
so  much  mass  that  the  liquid-gas  phase  transition  is  barely  driven  before  the 
point  of  maximum  current,  thus  allowing  very  little  time  after  the  phase  transi- 
tion for  heating  to  occur  and  the  temperature  to  rise.  The  highest  temperature 
is  not  the  ideal  case  for  laser  heating;  recall  that  the  inverse  bremsstrahlung 
absorption  coefficient  goes  as  Thus  increased  temperature  is  an  adverse 
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Figure  11.  Maximum  Temperature  Versus  Time  for  Aluminum  Exploding  Wires 
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effect  for  laser  heating  since  the  absorption  coefficient  goes  as  for 


less  than  the  critical  density.  Thus  it  is  important  to  hit  the  plasma  as  close 
to  the  critical  density  in  time  as  possible,  even  using  the  laser's  leading  pulse 
to  help  form  the  plasma  into  a good  heating  configuration.  It  is  convenient  in 
this  case  that  the  plasma  has  expanded  to  reach  its  critical  density  long  after 
the  maximum  temperature  is  reached;  at  this  point  the  plasma  has  cooled  off  and 
the  absor^'  .ion  coefficient  is  somewhat  increased.  The  expansion  velocity  for 
the  3-mil  wire  as  measured  from  the  plot  is  found  to  be  on  the  order  of  1.5 
centimeter/microsecond.  Recalling  that  the  density  goes  as  l/r^  cylindrical 
geometry),  the  3-mil  wire  would  have  to  expand  to  0.3-cm  radius  before  the 
density  is  at  the  critical  density.  This  corresponds  to  an  expansion  time  of 
0.2  microseconds  after  the  transition  which  allows  hydrodynamic  expansion  to 
occur.  This  seems  to  be  the  ideal  point  in  time  to  heat  with  the  laser.  For 
the  case  of  the  3-mil  wire,  which  seems  to  be  the  optimum  of  the  wire  considered, 
the  lapse  time  required  to  meet  optimum  conditions  is  0.8  microsecond  after 
initiation  of  the  capacitor  bank  which  explodes  the  wire. 


LASER  HEATING 


Given  the  parameters  of  the  preformed  plasma  as  indicated  in  the  initiation 
phase  of  this  study,  the  problem  now  becomes  what  heating  can  occur  from  a CO2 
laser  and  what  kind  of  scaling  exists  with  laser  pulse  length,  power,  and  total 
energy.  This  parameter  study  was  conducted  with  a one-dimensional,  laser  hydro- 
dynamics code,  DULAH2,  a modified  version  of  DULAH  (ref.  3)  as  described  in 
section  II  of  this  report. 


As  discussed  qualitatively  in  the  initiation  phase,  the  optimum  heating 
situation  exists  in  the  3-mil  wire  approximately  0.8  microsecond  after  bank 
initiation.  At  this  point  the  wire  is  a low- temperature  plasma  (0.5  to  1.0  ev) 


of  density  near  the  critical  density  for  10.6-ym  radiation  (N^  = 1.0  x 10^^). 


Beginning  with  an  aluminum  plasma  of  this  type,  DULAH  2 calculations  were 
made  with  a CO2  laser  pulse  like  that  obtainable  from  the  Air  Force  Weapons 
Laboratory  (AFWL)  cold-cathode  electron  beam-sustained  CO2  laser  (ref.  9).  This 
laser  typically  has  a gain-switched  pulse  of  50  to  80  nanoseconds  and  50  to  100 
joules  total  energy. 


A simulation  was  made  of  the  mode-locked  structure  (see  figure  12)  with  the 
same  total  laser  energy  in  a pulse  of  constant  power  level,  i.e.,  square  wave. 
The  electron  and  ion  temperature  profiles  for  these  two  simulations  are  compared 
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Figure  12.  Simulation  of  Mode  Li 
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in  figure  13.  Several  observations  can  be  made  from  this  comparison:  (1)  the 

maximum  electron  temperature  reached  in  each  case  is  only  slightly  different  and 
is  thus  more  a function  of  energy  and  power  than  pulse  shape;  (2)  the  maximum 
ion  temperature  reached  in  each  case  is  similar,  differences  are  due  to  ion  shock 
heating;  (3)  the  maximum  ion  and  electron  temperatures  occur  earlier  in  time  than 
the  cutoff  of  the  pulse  due  to  the  hydrodynamic  expansion  which  results  in  lower 
densities,  lower  densities  cause  lower  laser  absorption  in  the  electrons  and 
lower  energy  exchange  to  the  ions  from  the  electrons;  (4)  the  time  integral  of 
the  electron  temperatures  in  each  case  is  very  similar;  (5)  in  each  case  the  ions 
begin  cooling  off  before  the  electrons,  while  the  electrons  are  still  being 
heated  by  the  laser,  the  ions  are  cooling  because  of  the  hydrodynamic  expansion. 

In  conclusion,  the  mode-locked  pulse  does  not  yield  significantly  different 
results  from  the  "somewhat  ideal"  square-wave- type  pulse. 

The  next  question  to  arise  is  how  the  laser  interaction  is  affected  by  the 
scaling  of  laser  power.  Figures  14  through  16  show  some  comparisons  of  the  plasma 
dynamics  and  the  ion  and  electron  temperatures  for  laser  powers  of  three  different 
orders  of  magnitude.  These  laser  powers  are  in  units  of  ergs/cm/ rad/sec  which  is 
an  azimuthally  symmetric  laser  flux  incident  upon  a cylinder  in  the  r direction. 
For  the  plasma  of  interest  in  these  exploding  wire  calculations,  these  fluxes 
correspond  to  300  joules/cm,  3,000  joules/cm,  and  30,000  joules/cm  for  the  10^^, 
10^^,  and  10^®  cases,  respectively.  Figure  14  is  a comparison  of  the  plasma 
outer  radius  versus  time  for  these  three  cases.  This  shows  clearly  that  the 
effect  of  the  increased  power  is  to  blow  the  plasma  apart  on  a much  smaller  time 
scale.  The  slopes  of  these  curves  represent  the  velocities  of  the  outer  boundary 
of  the  plasma.  Estimates  of  these  velocities  for  the  various  cases  are  given  in 
table  2.  Figure  15  is  the  ion  temperature  versus  time  for  these  power  levels. 

This  shows  some  interesting  effects:  (1)  the  ion  temperature  peaks  earlier  in 


time  with  higher  laser  power;  (2)  the  width  of  the  temperature  versus  time  peak 
narrows  considerably  in  time  with  higher  power;  (3)  there  appears  to  be  some 
optimum  laser  power  for  heating  the  ions  because  the  maximum  ion  temperature 
reached  in  the  10^®  and  10^’’  cases  are  each  lower  than  the  lO^'^  case.  This  will 
be  borne  out  in  later  comparisons  as  well.  From  figure  16,  electron  temperature 
versus  time  for  these  power  levels,  the  following  conclusions  can  be  drawn:  (1) 

the  electron  temperature  appears  to  have  some  direct  scaling  relationship  with 
the  laser  power  (this  will  be  exploited  more  in  a later  comparison);  (2)  the  peak 
in  the  electron  temperature  appears  earlier  in  time  with  the  higher  laser  power. 
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Figure  14.  Outer  Plasma  Radius  Versus  Time 
as  a Function  of  Laser  Power 
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Figure  15.  Ion  Temperature  Versus  Time  as 
a Function  of  Laser  Power 
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Figure  16.  Electron  Temperature  Versus  Time 
as  a Function  of  Laser  Power 
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Table  2 

VELOCITY  OF  EXPANSION  OF  OUTER  PLASMA  BOUNDARY  FOR  VARIOUS 

VALUES  OF  LASER  POWER 


Laser  Power 
(ergs/cm/rad/sec) 

Plasma  Velocity 
(cm/sec) 

1016 

3.3  X 10^ 

IQi^ 

7.3  X 10^ 

IQiP 

1.3  X 108 

falling;  (3)  there  appears  to  be  little  difference  in  the  late-time  (beyond  75 
nanoseconds)  behavior  of  the  temperature  as  a function  of  laser  power. 

These  kinds  of  comparisons  raise  a question  which  is  the  most  pertinent 
question  which  can  be  asked  of  present  simulation  experiments;  that  is,  how  do 
these  quantities,  most  importantly  ion  and  electron  temperature,  scale  with  laser 
energy?  In  other  words,  are  these  experiments  worth  scaling  up  for  simulation 
purposes?  The  answering  of  this  question  can  be  helped  somewhat  with  figures  17 
and  18.  Figure  17  portrays  the  scaling  of  maximum  electron  and  ion  temperatures 
versus  laser  energy  for  two  different  laser  energy  for  two  different  laser  pulse 
lengths.  That  is,  holding  the  pulse  length  constant  and  increasing  the  energy 
in  the  pulse,  how  does  the  maximum  temperature  reached  scale  with  laser  energy? 
Figure  17  uses  the  laser  pulse  shape  shown  in  figure  12  with  pulse  lengths  of  10 
and  80  nanoseconds.  These  two  pulse  lengths  are  used  primarily  to  give  some 
qualitative  comparison  of  the  effect  of  changing  the  pulse  length.  However, 
numerical  relationships  in  the  linear  approximation  (of  a log-log  relationship) 
have  been  obtained  for  the  scaling  in  these  two  cases  for  the  electron  temperature 
alone.  The  scaling  laws  are  of  the  following  form: 

C*  = a(i\)'  (35) 

is  the  maximum  electron  temperature  in  electron  volts,  is  the 


where  T, 


•I  ^ 

total  energy  in  the  laser  pulse  in  joules,  and  A and  B are  empirical  constants 
given  in  table  3.  It  can  be  observed,  at  this  point  as  well,  that  the  relation- 
ship of  the  ion  temperatur'e  versus  laser  power  hinted  at  in  the  discussion  of 
figure  15  can  now  be  clearly  demonstrated.  There  is  a laser  power  or  total 
energy  beyond  which  the  maximum  ion  temperature  produced  declines  with  additional 
incident  energy.  The  hydrodynamic  expansion  occurs  on  shorter  time  scales  than 
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Figure  18.  Maximum  Electron  Temperature  Scaling  Versus  Pulse  Length  for 
Total  Laser  Energy  of  0.1,  1.0,  and  10  Kilojoules 
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Table  3 

PARAMETERS  FOR  MAXIMUM  ELECTRON  TEMPERATURE  VERSUS  LASER 
ENERGY  FROM  EMPIRICAL  FITS 


Laser  Pulse  Length 

10  nsec 
80  nsec 


0.375 

0.500 


the  heating  via  shock  hydrodynamics  or  ion-electron  collisional  exchange.  The 
plasma  is  effectively  blown  apart  as  far  as  the  ions  are  concerned  and  the  limit- 
ing ion  temperature  is  thus  reached  at  some  intermediate  laser  power. 

Figure  18  portrays  some  information  which  is  directly  related  to  that  in 
figure  17.  The  contours  plotted  in  this  case  are  contours  of  iso-energy  for  the 
laser.  This  shows  the  effect  that  shortening  the  laser  pulse  at  constant  energy 
has  on  the  maximum  electron  temperature  reached.  This  can  also  be  used  to  deter- 
mine whether  it  is  advantageous  to  shorten  the  laser  pulse,  given  a certain  energy 
change.  For  example,  shortening  the  laser  pulse  from  80  to  10  nanoseconds  with  an 
energy  reduction  from  1 kilojoule  to  100  joules  would  result  in  a slightly  higher 
maximum  electron  temperature,  whereas,  changing  from  an  80-nanosecond  to  a 30- 
nanosecond  pulse  with  the  same  energy  change  as  above  would  not  be  advantageous. 

3.  CONCLUSIONS 

The  most  important  results  of  this  research  are  twofold:  (1)  given  the  pulse 

length  and  total  energy  of  the  CO2  laser  used,  a good  estimate  of  the  plasma 
temperature  and  dynamics  can  be  obtained  within  the  liniitation  that  the  codes 
used  in  these  calculations  are  one-dimensional  and  ignore  some  important  two- 
dimensional  effects  due  primarily  to  the  extra  degree  of  freedom;  (2)  calcula- 
tions reveal  that,  at  least  in  the  40-  to  40,000-joule  range,  the  electron 
temperature  reached  in  the  plasma,  and  thus  the  radiation  output,  scale  at  least 
weakly  with  the  laser  energy  at  fixed  pulse  width.  The  amount  of  radiation  output 
and  the  conversion  efficiencies  are  not  taken  into  account  in  this  calculation  and 
are  important  quantities  which  remain  yet  to  be  determined  before  any  accurate 
assessment  of  the  simulation  possibilities  of  the  exploding  wire  can  be  made.  It 
would  appear  at  this  juncture  that  further  investigation  and  research  in  this  area 
are  indeed  justified  and  warranted. 
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SYMBOLS 


atomic  weight 

bremsstrahlung  radiation  term 
electron  thermal  conduction  term 
ionization  energy  term 

inverse  bremsstrahlung  absorption  coefficient 

electron  density,  cm"^ 

critical  electron  density,  cm"^ 

ion  density,  cm-^ 

electron  pressure 

ion  pressure 

ion-electron  collisional  exchange  term 
bremsstrahlung  radiation  term 
zone  radius,  cm 

thermonuclear  reaction  rate,  sec^ 
electron  temperature,  °K 
atomic  number 
average  ionization  levels 
speed  of  light,  cm/sec 
Boltzman  constant,  ev/°K 
electron  temperature,  ev 
electron  rest  mass 
electron  charge 
free  space  permittivity 
electron  internal  energy,  ergs/gm 
ion  internal  energy,  ergs/gm 
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Xg.i(Z) 


SYMBOL  (Cont'd) 


classical  Coulomb  logarithm 


laser  frequency,  sec 


plasma  frequency,  sec 


I 


ion-electron  collision  time 


laser  flux,  ergs/cm/radian/sec 


critical  optical  angle  of  incidence 


ionization  potential  to  remove  Zth  electron 


coefficient  of  electron  thermal  conductivity 
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